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Abstract 

^' 

CU The eigenvalue probability density functions of the classical random matrix ensembles have 

^ (— I a well known analogy with the one component log-gas at the special couplings /3 = 1,2 

and 4. It has been known for some time that there is an exactly solvable two-component 
log-potential plasma which interpolates between the (3 = 1 and 4 circular ensemble, and 
an exactly solvable two-component generalized plasma which interpolates between (3 = 2 
and 4 circular ensemble. We extend known exact results relating to the latter — for 
■ the free energy and one and two-point correlations — by giving the general {k\ + Ap- 

point correlation function in a Pfaffian form. Crucial to our working is an identity which 
expresses the Vandermonde determinant in terms of a Pfaffian. The exact evaluation of 
£^ ■ the general correlation is used to exhibit a perfect screening sum rule. 

j— i ■ 

1 Introduction 

1.1 The circular ensembles of Dyson 

In classical random matrix theory (see e.g. [6]), the probability density function (PDF) 

J_ TT | e **_ e <Y (1) 

for fj = 1, 2 and 4 plays a distinguishing role. Thus (JT|) for these values of /3 is the eigenvalue 
PDF for Dyson's circular ensembles of unitary random matrices {UnUJ^}, {Un}, {^iv^^} 
respectively. Here Un £ U(N), chosen with Haar measure and U D = Z2nU t Z^n, where 

Z 2 n = In® ^ q denotes the quaternion dual. The spectrum of matrices U 2 nU® n is doubly 

degenerate, which explains why (JT]) still involves iV eigenvalues for /3 = 4. 

Dyson [3] introduced the circular ensembles as alternatives to the Gaussian ensembles of 
Wigner. Unlike the latter, they are uniquely determined by an invariance property: their mea- 
sure is unchanged under left and right conjugation by real orthogonal, unitary, and symplectic 
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unitary matrices respectively. Significantly, it was found by explicit calculation that in the bulk 
scaling limit the two-point correlation functions for the two classes of ensembles are identical. 
Extrapolating this to the conjecture that for a given (3, all statistical properties of the bulk 
states of the circular and Gaussian ensembles are identical was important at the time, as it 
allowed calculations to be carried out on the computationally simpler circular ensembles, with 
the aim of answering questions relating to Gaussian Hermitian matrices. 
A celebrated example is the statement that [6j eq. 18.160] 



2p b 1 ulk (2k+l;2s) 



(2) 



where p^ nYk (k; s) denotes the PDF that a gap of size s between two eigenvalues contains k 
eigenvalues in the bulk state at coupling 0. This was deduced by deriving for the circular 
ensembles that [T5|, eq. (8.157)] 

PN A {n; (-9, 9)) = p 2N ,i(2n + 1; (-9, 9)). (3) 

Here the notation PN t p(k; —(9,9)) is analogous to p h ^ Xk {k; s) but with the extra subscript A" 
referring to the total number of eigenvalues. 

The proof of the conjecture that for (3 = 1, 2 and 4 the bulk states of the circular and 
Gaussian ensembles are the same became possible after another landmark paper of Dyson [I]. 
In this work a formalism was given which allowed not just for the calculation of the 2-point 
correlation functions, but the general /c-point correlations pa.y This formalism was applied to 
the circular ensembles. In the bulk scaling limit, particle density equal to unity, it gave 



for (3 = 2, and 



P(k){xi, 



p (fc )(xi, ...,x k )= det 



sin TC (Xa 



n(xj - xi) 



3,1=1, 



(4) 



qdet 



sin27r(xj — x{) 
27r(xj - xi) 
d sin27rp(£j — xi) 



dxj 2ir(xj 



Xl) 



2^ 

sin27r(xj — x{) 
2n(xj — Xi) 



smx 



x 



dx 



(5) 



-I o,i=i, 



for /3 = 4 (a formula analogous to (JSJ) was also given for /3 = 1, but as it plays no role in our 
subsequent discussion we refrain from writing it down). The 2k x 2k matrix, S say, in fl5]), has 
the property of being self dual: S D := Z2kS T = S. It follows that SZ^ is antisymmetric 
so its Pfaffian is well defined. The quaternion determinant in (jSJ) is then given in terms of 
a Pfaffian according to qdet S = Pf(SZ^). The method of [1] was subsequently applied to 
the Gaussian ensembles [H], allowing in particular the verification that the bulk statistical 
properties coincide with that of the circular ensemble. 



1.2 From the circular ensembles to two-component log-gases and 
trial wave functions 

As emphasized by Dyson (3], for general /3 > (pQ) has the interpretation as the Boltzmann 
factor for a classical statistical mechanical system with A^ particles on a circle interacting via 
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the pair potential — log |e* 9 — e % ^\ (one-component log-gas). Sutherland [18] showed that (JTJ) 
is the absolute value squared of the ground state wave function for the quantum many body 
system specified by the Hamiltonian 

2^ de yA\2 J 2s sin 2 7r(^-^)/2' {> 

3=1 3 l<j<k<N K °" 

These alternative interpretations suggest two-component generalizations of (JTJ which interpo- 
late between the classical couplings (3 = 1, 2 and 4. 
One such generalization is (5] 

Ni N 2 

1 [ |e^-e^'| If le^-e^^JJJJle^-e^l 2 . (7) 

l<j<k<Ni l<a<(3<N 2 j=l a=l 

It is the Boltzmann factor for a two-component log-potential plasma system on the circle with 
N\ particles of charge +1, and N 2 particles of charge +2. It is also the ground state wave 
function for a two-component generalization of ([6]) [llj. For N 2 = (Ni = 0) (JTJ) reduces to (111) 
with = 1 (0 = 4). Moreover, the general (ki, fc2)-po m t correlation function can be computed 
in a quaternion determinant form [SJ eq. (6.168)] which interpolates between the quaternion 
determinant forms for (3 = 1 and (3 = 4. And in a very recent development the Gaussian 
analogue of (JTJ) has shown itself to be exactly solvable [16]. 
A second such generalization is [TJ 

_^ ]"J | e ^- e ^| 2 ]"] | e ^-e^| 4 ]Jf||e^-e^| 2 . (8) 

2 l<j<k<Ni l<a<P<N 2 j=la=l 

This describes a Boltzmann factor of a so-called generalized plasma: the couplings (exponents 
on the products of differences) g^R, god 9rg f° r the two species (R)oman and (G)reek no longer 
satisfy gnngcc = 9rg as they do in (JTJ). With the complex exponentials replaced by general 
complex coordinates, this is of interest as a trial wave function in the anomolous quantum Hall 
effect [S] , and without such replacement as the exact wave function for a certain supersymmetric 
tJ system in one-dimension [QQ1]- For N 2 = (JVi = 0), flSJ) reduces to CQ) with (3 = 2((3 = A). 
However, unlike the situation with the two-component system (JTJ), there is no existing formula 
in the literature for the general (k\, /c2)-point correlation function, nor are there any exact 
calculations carried out for the Gaussian analogue of (JHJ). 

Two exact solvability features of (JHJ) are in the existing literature. One is the exact evaluation 
of the normalization, 

and the other is an exact evaluation of the two-point correlations. The latter, in the bulk 
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scaling limit with density of Roman species pr and Greek species p G , read j7] 



2 sin 2 Tr/jjja; 
Phr(z, 0) = p H - 



(7nc) 2 

-PGPR I dt I ds ( e 2-Ki PR x{s-l)-2*i PG xt _ e 2m PR xs+2^ PGX t\ ^ 

Jo Jo 2p G t/p R + l\ J 

+ (2t + P Jla)(L r IPo) C ° S ^ PGX(t + ' + PRlPc) )) (U) 

Pi? G (x, 0) = p R p G - 2p/jp G / dt I ds ^ ct + PrS cos (2np G xt - 27rp jR x(s - 1)) , (12) 

io Jo 2p G t + p R V / 

which can of course be equally as well written as products of one-dimensional integrals. In 
the case (pc,pn) = (1,0) the first of these reduces to @ with k = 2, while in the case 
(Pg,Pr) = (0, 1) the second reduces to (jSJ) with k = 2. 



1.3 Our aim 



In light of the above discussions, two challenges present themselves. 

1. To compute the (k\, k 2 )-pomt correlation function corresponding to the generalized two- 
component plasma (|S]) in a form which relates to (pEJ) and (151) . 

2. To similarly compute the (ki, A^-point correlation function for the Gaussian analogue of 

We will see how application of a Pfaffian structure only very recently introduced into random 
matrix theory [T7j allows us to solve the first of these problems. In the Appendix we generalize 
this Pfaffian formula to a form suitable for application to the second problem, but we leave the 
required subsequent working to a later work. 

The previously computed two-point correlations f|T0|) -f fT2|) hold in the bulk scaling limit. We 
obtain the corresponding (k\, /c2)-point correlation functions by taking the bulk scaling limit of 
our finite system results. In the work [7] some special screening properties of the bulk scaled 
truncated two-point functions, 

/oo pec /*oo 

Pr R (x, 0) dx = -p R , / Pr G (x, 0) dx = 0, / p GG (x, 0) dx = -p G (13) 
-oo J~ oo J— oo 

were noted. Knowledge of the {k\, A^-point correlation function allows these sum rules to be 
generalized so as to relate to the truncated form of the general correlation functions. Further- 
more, we give a discussion relating to the correlations in the large pr or large p G limit. 
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2 Correlations for the circular generalized plasma 
2.1 Generalized partition function 

The method used in [7j to deduce the exact two-point correlations (|T0l)-(IT2l) was to write ( 
as the product of a Vandermonde and confluent Vandermonde determinant, by noting that 



n (** ~ z ^ = det ^j i ]i,fe=i,-,iv 1 

l<j<k<Ni 



H (z k - z^** = det 
i<j<k<N 1 +N 2 



J fe=l,...,JV,+2JV 2 



ik - \)z)-\ 



j=Af 1 + l,...,JV 1 +JV 2 
fe=l,. ..,N 1 +2N 2 



(14) 



(15) 



where = 1 (j = 1, . . . , iVi), ^ = 2 (j = iV 2 + 1, . . . , iVi + iV 2 ), respectively. 

The most crucial strategic step of the present study is to replace f }T4|) by the Pfaflian identity 



l<j<k<Ni 



[Z 



Ni/2 m/2s 2 



J 



Zj~ Zj 



j,k=l,...,Ni 



(16) 



We give a new proof of this identity, and a generalization which is of potential future use in 
random matrix theory in the Appendix. This structure was introduced into random matrix 
theory in the recent work of one of us [IT]. It allows us to express the generalized partition 
function 



Z Nl ,N 2 [u,v] := 

* n 

i<j<k<Nt 



1 



C7(JV!,JV 2 ) J 



2tt 



d9 1 u(6 



2n p2tt 

d6 Nl u(0 Nl ) j #iu(0i) 



e ie k _ e iej\2 



n 

l<a</3<N 2 



'0 JO 

Ni N 2 

|g#a _ e i<f>/3 |4 TJ TJ | e «0j _ g^a |2 

i=l a=l 



2tt 



C?0JV 2 W(0AT 2 ) 



(17) 



in terms of a Pfaflian. 

Proposition 2.1 VKzt/i 2; := e l6 ' , and {Pj-i{z)}j=i,2,... an arbitrary set of raonic polynomials, 
Pj-i(z) of degree j — I, let 



2- 



2tt 



a^j = / / d9 2 u(9 2 )z^z^ N AZ2 - x \ '- Pj-i{zi)p k -i(z2) (18) 



-JVi/2 _ -iVi/2^2 
fl_ 

z 2 — 2^ 



2tt 



„(fl) Z -(«+^( ft „ 1 (,),/ t _ 1 ( z ) - ^(^(2)) (20. 



(19) 



In terms of this notation, and with [C fc ]/(C) denoting the coefficient of ( k in the power series 
expansion of /(C); for Ni even we have 



Z Nl ,N 2 [u,v} 



Nx\N 2 \ 
C(N U N 2 



iC^PnCaM + hMkk^m- 



-2N 2 - 



(20) 
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Proof. Using the identity 



\ e i6 k _ e i6 j \ 



-ie 



-i(6j+e k )( je k _ in- 



valid for < 9j < 9 k < 2n the absolute values appearing in the integrand can be eliminated, 
and we obtain 



Zni,n 2 K v] :-- 



2t 



dOiuiOi) 



2t 



2~ 



hv((pi ■■■ 



2tt 



Ni 



xn^n^t +2Ni - 2) n fe 1 -^ 1 ) n 

1=1 a=l l<j<k<N 1 l<j<k<N 1 +N 2 

Consider the Pfaffian formula ffTB]) with the replacements (zj t— > zj l ). It can rewritten 



d(j) N2 v((f) N2 ) 
(21) 



n (^-^) 

l<j<k<Ni 



1 



2^/2(^/2)! 



Asym | | 



iVi/2 / -JVi/2 _ -JVi/2x 2 
Z 2j Z 2j-1 J 



Z _1 - 

Z 2j Z 2j-1 



(22) 



where Asym denotes the operation of antisymmetrization in {zi, . . . , zn x }- Substituting ( 122]) in 
(T2~T1) we see from the fact that the final product of differences is antisymmetric in {z±, . . . , zn ± } 
that all terms in the Asym operation contribute equally. Thus 



Z Nl ,N 2 [u,v} 



C(N U N 2 ) J 



2tt 



A'2 



X 



-(jVi+2JV 2 -2) 



d9i u( 



Ni/2 f -Ni/ 2 



2- 



2- 



dO Nl u(6 Nl ) / d(j>x v(4>i) ■ 



2- 



1=1 



a=l 



n 



2.7 



— Z, 



-iVi/2\2 
2j-l ) 



Z 2j Z 2j-1 



n ^ - z i- 

l<j<k<Ni+N 2 



d(j) N2 v(4>N 2 ) 
(23) 



where C{N X ,N 2 ) = 2 N '/ 2 (N 1 /2)\C(N 1 , N 2 )/N 1 \. 

Next, we observe (as is standard) that for {pj(x)}j=o,i,2,- a se ^ °f mon i c polynomials, Pj(x) 
of degree j, the determinant in (j2J) can be rewitten to read 



det 



\Pk-x(zj)] 3=1,...,^ 

fc = l,...,iV 1 +2JV2 

Pfc-l(fj) 

(k ~ l)p' k -i(zj) 



j=N 1 +l,...,N 1 +N 2 
k=l,...,N 1 +2N 2 . 



N, 



N ± +N 2 



Yi £ (q)Y[pq(i)-i( z j) n p<3(j)-i( 2 j)PQ(f)-i( z j). 

QG5JVJ+2JV2 



z=l 



J=JVi+l 



where the equality follows from the definition of a determinant as a sum over permutations. 
Substituting this for the final term in (12"3"|) . we see that the integrations factorize down to one 
and two dimensional integrals, and we obtain 



1 



C{Nx,N 2 

Ni\N 2 \ 
C(N U N 2 



Ni/2 N 1 /2+N 2 

£ (q) n a Q(«-i),Q(2oM n 6q(«-i),o(2I)H- 

QeSiv 1+ 2iv 2 J=l !=AT 1 /2+l 

JVi/2 Vi/2+V 2 

*e(<5) II °P(a-i),P(2oM II 6 Q(2/-i),Q(2oN- ( 24 ) 
QeS Nl+2 N 2 j=l i=Ni/2+i 
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Here 

Jo 

and the asterisk denotes that the permutations are required to satisfy P(2l) > P(2l — 1) for 
eachZ = 1, . . . , N t /2 + N 2 and P(l) < P(3) < •■■ < P(N 1 + 2N 2 -1) (imposing these conditions 
is how the second equality is obtained from the first). But by writing out the Pfaffian in fl20]) 
as a sum over permutations according to the general formula 

TV 

pfc= Y.*<Q)Ii 

CQ(2l-l),Q(2l) 

Q £ 5 2JV 1=1 

for C = [cjk]j,k=i,...,2N antisymmetric, we see it reduces to (|24|) . □ 



2.2 Skew orthogonal polynomials 

Let R denote a particular permutation of {1, 2, . . . , iVi + 2iV 2 }. Let us rearrange the rows and 
columns of the Pfaffian in (1201) to have the order of R, and thus 



Z N1)N2 [U,V) = J^j^s [C Nl/2 ]Pf [CaR{j),R{k)[ U ] + b R (j),R(k)Mj,k=l,...,N 1+ 2N 2 - (25) 

We seek a choice of the permutation R, and a choice of the polynomials {pj{z)}, such that the 
skew inner product C a R(j),R( k ) + ^R(i)M k ) nas ^ ne skew orthogonality property 

tn rn _i_ /i hi / r i-i> (^,"i) = (2j - l,2j) 

C^),*M [1] + &k ( „),*m [1] = | Q) otherwis6) (26) 

for P(n) < P(m). 

Proposition 2.2 For j = 1, . . . ,N±/2 set 

R(2j -1) = N 2 + j, R(2j) = N 2 + N 1 -j + l (27) 
and for j = N x /2 + 1, . . . , iVi/2 + N 2 set 

R{2j - 1) = j - iV x /2, P(2j) = 2iV 2 + 3iVi/2 - j + 1. (28) 

C/ioose 

Pi(z) = zi {j = 0, 1, . . . , iVi + 2iV 2 - 1). (29) 
VFe t/ien /iai>e that the skew orthogonality property &b)) holds with 

f (2tt) 2 C + 2 7 r(iV 1 + 1 - 2j), j = 1, . . . , iVa/2 
^ \ 2^(2^ + 2^! + 1 - 2j), j = N 1 /2 + l,...,N 1 /2 + N 2 . [ > 
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Proof. With w = e*^, let us write 



V27T ( z -Nl/2 _ -JVi/2\2 

*i(*)'-= ™~ N A y -r (31) 
Jo z a? 



Use of the expansion 



2 1 — W 1 



iVi/2-1 

^ (j*- 1 )^" 1 )^/ 2 - 1 -* (32) 

p=0 



shows that 



$ f 27rsgn(iV 2 + iV 1 /2-l/2-j>^ 1 -^-^, j = N 2 , . . . , iV 2 + JV X - 1 (33) 

5 0, otherwise. 

But according to the definition (fTSl) 

a i|Jfc [i]= r z-^-^-^de. 

Jo 

Making use of (|33|) we therefore have 

a- [11 = / ( 27r ) 2s S n ( Ar 2 + iVi/2 + 1/2 - .7)^-1-^-2*3,0, j = N 2 + 1, . . . , N 2 + 
J,fc ^ J \ 0, otherwise. ^ ' 

On the other hand, from the definition (j!9p and (1291) we see that 

6 i>fc [l] = (fc-j) / 2 %- (Ari+2Af2)+J+fc - 1 ^ = 27r(A:-j)Wi-^-2* 2 ,o. (35) 

In light of fl3^|) and fl35|) . we see that the skew inner product in ([26]) is nonzero for i?(j) + 
R(k) = Nx + 2N 2 + 1 only. But ([27]) and ([28]) tell us that for j < k, this equation holds if and 
only if j = 2j' — 1 and k = 2j'. In these latter cases we have 

C a R(2j-l),R(2j) + bR( 2 j-l),R( 2 j) 

= ((27r) 2 Csgn(iV 2 + N,/2 + 1/2 - R(2j - 1)) + 2tt(2N 2 +N x + 1- 2R{2] - 1))) 

for j = l,...,N x /2 and 

Ca R ( 2j -i), R ( 2j )[l} + b R{2j _ 1)M2j) [l] = 2vr(2iV 2 + N x + 1 - 2i?(2j - 1)) 

for j = JVi/2 + 1, . . . , JVi/2 + N 2 , which establishes the result. □ 
According to (121)]) 

AM AM N2/2+N2 

z ™M=mk>^ n ^- 

Substituting ( 130]) and simplifying we see that for Zn 1} n 2 [1, 1] = 1 the normalization must be 
given by <jHJ). 
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2.3 PfafRan form for the correlations 

Consider a two-component system of (R)oman and (G)reek particles. Let the one-body po- 
tentials u and v be associated to each species respectively, and suppose the corresponding 
generalized partition function has the form 

Z N [u, v] = Pi[aj )k [u] + bj,k[v]]j,k=i,...,N- (36) 

Here N is assumed even and aj )k [u\, bj tk [v] are antisymmetric in j, k. 

Generally for a two-component system the (ki, k 2 ) -point correlation function can be com- 
puted from the generalized partition function by functional differentiation according to the 
formula 



P(fci,fc 2 )( x l; • • • > x ki) Vl, ■ ■ • , Uk 2 ) 

I fiki+k-2 
= — —. r —. r-z—. r —. -Zjsr[u,v] 

Z N [l,l]ou(x 1 )---du(x kl )dv(yi)---dv{y k2 ) 



(37) 

u=v=l 



Our first point is that with some additional structure, the generalized partition function (I36I) 
substituted into (1371) yields a Pfaffian of size 2(ki + k 2 )- 
Thus we suppose that 



a jA u } = 2 J u ( x )^j R:) £r( u 4>1 R ^)[ x } dx 

b j>k [v) = 2 j v{y)(4,f\y)e G ^f%]-4>f\y)e G {4,f%})dy, (38) 



where, with s G {R, G}, the e s are linear operators. We see that ( j2Qj) is of the form ( l38l) with 
I = V = [0, 2vr] and 

^f\ x ) = z - N > PRa+iyi (z), z = e 2 ™ (39) 

e R (fM = i / r ^f(z 2 ) dz 2 (40) 

2 Jo z 2 - z- 1 

i>f\x) = z-^ +2N >- 2 »*p Rij+iyi (z) (41) 

^ = ItM^ (42) 

Furthermore we know from Proposition 12.21 that there is a choice of {pj{x)} which skew diag- 
onalizes the inner product. Let us assume that the ijy?' are chosen so that the inner product 
implied by (I3"6"|) with u = v = 1 is similarly skew-diagonalized. In this setting the explicit form 
of the corresponding (ki, A; 2 )-point correlation function (13"7|) has been computed by Mays 
(see also [2]). 
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To present the result, let 

AT/2-1 



N/2-1 

D SUS2 (n, v ) = 2 -\4f(^4^i(v)-^UMf(r]) 

j=o j 

N/2-1 

i SUS2 M = 2 Yl -(tMfK^sM-iiM - e s A4^i)(^sA4f](v) 



3=0 ] 



+ < 2 z- 1 - w- 1 



si = s 2 = R, 
otherwise. 



Also, set 



~S JuS2 (n,r}) -D SljS2 (n,r])~ 
(Note that K GR (y,x) = —K RG (x,y).) In terms of this notation, we have from [13] that 



(43) 



(44) 



1 



P{k\,k 2 )\ X li •••) x ki'i 2/1) ••■) ) — ~ r v y, 2--\- 



Nl/2+N 2 



xPf 



) W 
K G R(yi,Xj) K GG (yi,y, 



Zj 2(k 1 +k 2 ) 



i,j=l,...,k 1 
I ,m— 1, . . . ,/c2 



(45) 



An important point is that the expression in the numerator of ( 1451) is a polynomial in (. 
Expanding the Pfaffian for large ( will therefore allow us to determine what powers of ( in the 
expansion of rTil 1 / 24 ^ 2 r '-i wm combine to contribute to [C^ 2 ] ■ 



Proposition 2.3 We have 



lim Pf 



Pf 



KRR^Xi.Xj) K RG (xi,y n 
K GR (yi,Xj) K GG {y h y v 

K RR (xi,Xj) K RG (xi,y m ) 
K GR (yi,Xj) K GG (yi 

) ym) 



7- 1 

Zj 2(k 1 +k 2 ) 



i,3=l,-..,fcl 
Z,m — 1, . . . ,k 2 



7- 1 

Zj 2(k 1 +k 2 ) 



i,j=l,...,hi 

I ,m — 1 , . . . ,k 2 



where 



K S1 ,s 2 (^,V) ■= j im K s 1 ,s 2 (^,V) 

C—>oo 



(46) 



(47) 



is independent of(. Consequently 



P(ki,k 2 )(xi, ■ ■ -,x kl ;yi, ■ . -,yk 2 ) 



Pf 



K-RR^Xi^Xj) K RG (x i} y, 
K GR (y h Xj) K GG (y h y 



rn ) 
rn ) 



Zj 2(k 1 +k 2 ) 



i,j = l,...,k 1 
I ,m— 1, . . . ,k 2 



(48) 
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Proof. Assuming (j46p . the Pfaffian can be expanded in inverse powers of (, with leading term 
a constant. But from ( 130]) Yl^{ 2+N2 H-i is a polynomial of degree ( Nl / 2 and so the operation 
[C^ 2 ] i n the numerator of (|44|) is equal to ( 146|) times the denominator, and (148|) follows. It 
remains then to verify (jlSD- 

In fact substituting f[311 - (fl2"l) in (H3l it is immediate that all the leading large ( forms are 
independent of ( except for Irr(x,i/). Noting from fl27|) . fl33l) . fl39l) and f|40|) that 

(*)i^_J -K^ 1 ^ 1 , j = 0, . . . , AT a /2 - 1 



0, j = iVi/2, . . . , JVi + N 2 - 1 



eM (R) ](x) _) j=0,...,N l /2-l 



we have 



l R , R (x,y) = -- g (27r)2c + 27r(iVl + 1 _ 2j) ^ l * ^ 3 ) + 2 • 

This would appear to be proportional to ( for £ large. But it is easy to verify that in fact the 
proportionality constant vanishes, and so 

This like all the other entries of the Pfaffian is independent of ( as ( — » oo. □ 



2.4 Bulk scaling limit 

The probability density function (jSj) is defined on a circle. By scaling the angles 9j i— >■ 2nXj/L, 
4> a i — ^ 2iry a /L the system can instead be interpreted as defined on an interval of length L 
with period boundary conditions. The thermodynamic limit is then specified as L — )■ oo with 
N\/L =: pr and N 2 /L =: pc fixed. Applied to the correlation functions it gives the bulk scaling 
limit, which from [7] has the explicit forms (fT2l) in the two-point case. 

Using Proposition 12.31 it is a straightforward exercise to compute the bulk scaling limit for 
the general {k\ + fc 2 )-point correlation function. These involve the quantities 



^ ulk (a;, y)=p R f e 27ripR{x - y)t dt 
Jo 

r l/2 

l£ ulk (x, y) = p\ \ te 27TipR{x ~ y)t dt 

' J -1/2 

?bulk/ „\ _ PG / ( -2-Ki(p G t+ PR /2)(x-y) 







S G uik (x, y) = Y I ( e- Zm{pGt+pR/Z)(x - y> + c.c. ) dt 



I^x, y) = Yj (pct + P/?/2) [e- 2m( - PGt+pR/2)( - x - y) - c.c.J dt 
D hulk (x, y) = p G j - t ^ - (e-iMPGt+Pn/Wx-v) _ cx .) dt ( 51 ) 
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Proposition 2.4 Let 

p\k^k 2 )i x i,---, x k 1 ;yi,---,yk 2 ) 



We have 



^lirr^ J p (felifc2 )(2vrxi/L, . . . , 2irx kl /L; 2iryi/L, 2ny k JV) 



N 1 /L=p R 
N 2 / L =PG 



P(ki,k 2 )( Xl 
= Pf 



) x k\ i V\i 
-bulk 



,Vk 2/ 

bulk/ 
RG \ 



T^bulk/™ ™ \ 7>-bulk/™ . 
^ RR V d >l-> d >3) ^RG rn S/n 



-1 

2(fei+fc 2 ) 



(52) 



where 



with 



bulk (//, 77)' 



obulk 



s£ ulk (^): 



-^ k (*,y) 



J^ k (x,y) = -e-^ (x+ ^I h R ulk (x,y), 
-e nipliX D hulk (x,y), 



RR 

-D 



-e wipR(x+y) D bulk {x,y) 
S R t(x,y) = e m ^S G (x,y) 



bulk/ 
RG \ 



fbulk 
l RG 



(x,y) = -e-^I b R * lk (x,y) 



S GR (y,x) = e -^ i PR x e lxipRV S 

jbulk 



I^(x,y) = 

'GO ( X ' 



^bulk 



R{x,y), 

RG (y,x), S G u G \x,y) 
-D^ k (x,y), ± gg 



-D 



gr i x -> y) - 

S G ulk (x,y), 



DT(y,x) 



I^ k (x,y) = -I G ^ k (x,y). 



(53) 



Proof. We see from (125)1 . (150)1 . |j59l) (PT2D. together with flUD that the summations in (03 
range over j = 0, 1, ... , N\/2 in the cases that the summands involve e R \ipj ], and over j = 
Ni/2, Ni/2 + 1, . . . , N\/2 + N2 — 1 otherwise. In the former range 



^if(x) = z\ 



e G l4?](x) = -(--f + l+j); 
r J = (2vr) 2 C + 2vr(iV 1 + l-2j), 
while in the latter range 



-Ni/2- j 



x 



x) 



z 



-Ni+j-Ni/2 



z- 



J-JV1-JV2+I 



e G [^f](x) = -(j-iV 1 -iV 2 )^ 
r i = 27r(2iV 2 + 2JV 1 -l-2j), 



j-Ni-N 2 



4°Ux) 



z 



Ni/2-j 



(«) 

2j+l 



;r 



2 



JVi-l-j 



E-R^i+l]^) = 3 
= i(^l_ <7 -)^i/2-i-l 

2 V 2 ; 



(54) 



2j+l [ 



Z 



N 2 +3N 1 /2-j-l 



..(G) 
; 2j+l» 



Ni_+N 2 -j 



e G [^g 1 ](x) = -(iV 1 + iV 2 - J > 



2; 



ATi+JV2-i-l 



(55) 
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In relation to the cases that the summand contains e R [ipj ], appropriate use of ( |54l) sub- 
stituted in the first and third formulas of (|43|) . or (I50I) . shows that the summations are in fact 
Riemann sums. Simple manipulations then give 

2-7T f 1 

lim —S RR (27rx/L,27ry/L) = p R / e 2 ^ ^ dt 

L^oo L Jq 

lim -^-S gr (2itx/L, 27ry/L) = p R e^ ipRy / e -^ ip ^ x ' y)t dt 

L^OQ L Jq 

lim ^I RG (27rx/L,27ry/L) = \p\e~ mpRX te 2mp ^ x ~ y)t dt 

L-^oo I J -1/2 

2tt r l/2 

\im—I RR (2irx/L,2%y/L) = -p 2 R e-™ pR( - x+y U te 27Tip ^ x ~ y)t dt (56) 

L^oo L 2 J_ l/2 

(Ft) 

Similarly, in the cases that the summand does not contain e R [ipj ] we find 
lim ^S GG (2nx/L,2ny/L) = ^ / f e -^C«?*WO0»-v) + CiC .) ^ 
lim ^S RG (27ix/L,2ny/L) = ^- e nipRX [ ( e^t+^M + c . c dt 

L— >00 Z/ 2 Jq \ / 

hm ^/ GG (2vra;/L, 27ry/L) = Ap G t + Pfl /2) ( e -^*W9<*-i0 - c .c.) rft 

L->oo L z 4 Jo V / 

lim 2nD GG (2itx/L, 2ny/L) = p G [ ( e -^(PGt+ PR /2)(x- y ) _ \ df 

l->oo J p G t + p R /2\ J 

lim 2tiD rg (27cx/L, 2WL) = p G e"^ x / f ^(Pof+P^M _ c c \ ^ 

:-+oo J p G t + p R /2\ J 



lim 27iD RR (2nx/L, 2ny/L) = p G e mp ^ x+y ^ [ ( e ^{PGt+ PR /2)(x~ y ) _ \ df 

l^oo J p G t + p R /2\ J 

(57) 

Upon recalling the definitions (15"TI) . and making use of the fact that the Pfaffian and thus 
correlation function is unchanged if the element D S1S2 in (jUJ) is multiplied by the scalar L say, 
and I Sl . S2 i s multiplied by 1/L before the limit L — > oo is taken, we see that the stated forms 
follow. □ 

3 Properties of the correlations 

3.1 Explicit form of the two-point functions 

Denoting the bulk two-point function by p SlS2 (x,y) we see from fl52|) and fl53l) that 

Ps lS2 (x,y) = p Sl p S2 ~ (s si ,82(x,y)S S2>ai (y,x) + D ai>aa (x, y)I 8l>32 (x, y)\ 

Substituting the explicit forms of the quantities on the right hand side from Proposition 12.41 
we can check, after some minor manipulation, that there is agreement with the forms (fl2j) . 
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Furthermore, in the case pa = 0, we have from fl5Tl) that /)j ulk (:E, y) = 0. The structure of 
(1521 in the case k 2 = 0, and in particular the fact already alluded to below fl57|) that in the 
expansion of the Pfaffian each factor of / bulk is multiplied by _D buIk tells us that we can effectively 
set D huYk (x, y) = 0. The Pfaffian then reduces to a determinant, and this upon recalling the 
explicit form of S R nYk (x,y) from (15T1) is seen to be precisely the determinant formula (j4]) for 

n bulk 

P(k,oy 

We can check too that with Pr = the Pfaffian formula §5§ for is reclaimed. 



3.2 Sum rules 

Let us now turn our attention to the screening sum rules (fT3|) . In physical terms these say that 
upon fixing a Roman species particle, the other Roman species particles collectively redistribute 
to exactly cancel out the density of this particle. Similarly for the Greek species. This behaviour 
is to be contrasted to what happens in the two-component plasma specified by (JTj). There the 
analogue of (TL3I) reads 

/ (p+i,+i( a; >0) + 2p+i, + 2(^,0)J dx = -p +1 , / (p+2,+i(^0)+2p+ 2i+2 (x,0)J dx = -2p +2 . 

J —oo J— oo 

(58) 

In particular, unlike with (fT5|) . there is no special form for the individual integrals. In words 
the sum rules (|58p say that the effect of fixing a single +1 charge species is for both the +1 
and +2 charges to collectively respond to exactly cancel out this charge density, and similarly 
the effect of fixing a +2 charge species. 

Notice that (|T3|) and (1581) involve the truncated two particle correlations. In general the 
(ki, k 2 )-pomt correlation p^kl\ 2 ) ^ oes no ^ decay for large distances in any of its arguments and 
consequently cannot be integrated over R with respect to any one of its arguments. However, 
by adding suitable linear combinations of lower order correlations to P^ k k2 ) a quantity Po^J) 
- the fully truncated (hi, /^-point correlation - can be obtained, which has the property 
of decaying when any one particle coordinate is taken to infinity. The simplest case is when 
k\ + k 2 = 2 and we have for example p^ ^(x, 0) = p(i i i)(x, 0) — P(i,o)P(o,i)- A significant feature 
of the Pfaffian form (HH!) for the correlation functions is that it allows for p^ 1 ^ in the general 
case to be written in the structured form [6, Prop. 5.1.2] 

P(Xft5(*i, • • -> x *'>yi> ■ ■ ■ >v*) = (-i)* 1 ^ 1 E 



cycles 
length &i +&2 



/ \(o) 

X ^(^^^^(yfi.yia)^.^),^)^,^) ■ • " K »{*i kl+ik ys^iViH+i^Vii) ) ( 59 ) 

where {zi}i=x,... ) k 1 +k 3 = {^i}i=i,...,fci ^{yj}j=i,...,ki, s ( z ) refers to the species of coordinate z (thus 
s(x) = R, s(y) = G) and the operator ( )(°) refers to |Tr. 

For the two-component plasma the general truncated (k\ + A;2)-point correlation satisfies a 
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generalization of the sum rules fl58l) [6j (14.20)], 

P(fci+i,fc 2 )( a ' 1 ' • • • ' x kn x ] 2/i 5 • • • ) 2/fc 2 ) 

/oo 
pjk u k 2 +i)( x u ■■■, x k 1 ;yi,---,yk 2 ,y)dy 
-oo 

= -(fei + 2k 2 )pJ ki+1M) (xi, ...,x kl ,x;yi,..., y k2 ). (60) 

Here the factors of 2 are due to the charge density being the fundamental quantity. This sum 
rule is equivalent to the requirement that upon fixing +1 charges at xi, . . . , x kl and +2 charges 
at yi, ... ,yk 2 , the charge density of the system responds to exactly cancel these fixed charges. 

We know that for the generalized plasma, in the case of the two-point functions, the stronger 
sum rules (lT3"j) hold true. We will see that the existence of stronger screening rules is restricted 
to hi — 0, while more generally 



oo 

T 



P(jfci+l,fc 2 )( X l> • • • ' X kii x 'i 2/1' • • • > Vkz) dx 

/oo 
P(ki,k2+i)( x ^ ■■■> x k 1 ;yi,---,yk 2 ,y)dy 
-oo 

= -(fci + k 2 )pj kl+1>k2} (xi, . . . , x kl , x; j/i, . . . , y k2 ). (61) 

In contrast to fl60|) there are no factors of 2 indicating that screening is due to the total 
particle density rather than total charge density. The sum rule f loTj) is corollary of the following 
integration formulas for the matrix kernels (cf. [U Prop. 4.7]). 

Proposition 3.1 We have 

I (k rr (x 1 ,u)K rr (u,x 2 ) + K RG (x 1 ,u)K GR (u,x 2 )j du = K RR (x 1 ,x 2 ) 

J — oo 

K GG {yi, u)K GG (u, y 2 ) du = K GG (y±, y 2 ) 



K GR (yi, u)K RG (u, y 2 ) du = 

) 

/ (K GR (y, u)K RR (u, x) + K GG (y, u)K GR (u, x)jdu = K GR (y, x) 
/ (K RR (x,u)K RG (u,y) + K RG (x,u)K GG (u,y)j du = K RG (x,y). (62) 
Proof. The main formula required in the derivations is the generalized integral 

^ iks dk = 5{s) 

where 5(s) is the Dirac delta function. For example, using this formula we obtain 

K RR (x, y)K RR (y, x) dy = 



S R (x,x') 
-e- nip ^ x+x 'h R (x,x') S R (x',x) 



(63) 
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and 

K RG (x,y)K GR (y,x')dy = 2 



(64) 



Adding ( 1631) to ( 164|) and recalling ( 1571) we obtain the first of the sum rules in ( |62|) . 

The derivation of the remaining formulas in ( 1621) proceeds similarly. □ 

We see that indeed that after substituting ( 1591) in the LHS of (1601) and using the integration 
formulas ( 1621) that the RHS of (16T1) results. In the case of species G only, we similarly deduce 
the stronger result 



oo 

pf 0M+ i)(yu • • • , Vk 2 , y) dy = -k 2 pJo,k2)(yii Vk 2 )- (65) 



oc 



On the other hand, in the case of species R only, while p r RR satisfies the stronger sum rule in 
(IT3"|) . it follows from ( 16 3 p that we do not have for example that 



oo 

T 



is equal to —p]i R {xi, x 2 ). 



3.3 The large density limit of one of the species 

A possible limit of the two-component system is to take the density of one of the species to 
infinity, while keeping the density of the other species fixed. It turns out that this limit applied 
to the two-point correlations reveals a dramatic difference between the true log-potential plasma 
system (JTj) and the generalized plasma (JE]). 

For the generalized plasma, we see from ([TO]) and (|TT|) that 

hm o T (x 0) - Sin2w 

~Pr f dt [ ds - — ^ S ( e 2lTi PR x< y s ~ 1 )~' 2 ' ni PR xt — e 2irip R xs+2Trip R xt\ 

Jo J o 2t + 1 V / 
hm o T (x 0) - Sin27rpGX 

On the other hand, the explicit form of the two-point correlations for the true log-potential 
plasma system (j7j) [5J allows us to compute that 

lim p+ 1+1 (x,0) = 0, lim p+ 2 +2 (x, 0) = 0, 

where it is assumed a; / 0. Thus in the latter case the fixed density species become uncorre- 
lated when placed in a sea of the other species. This contrasts to the behaviour exhibited by 
the generalized plasma, in which for this setting the fixed density species still exhibits slowly 
decaying correlations. 
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Some insight into the reason for this is that, in the notation for the couplings of the gen- 
eralized plasma introduced below OH]), arguments based on the direct correlation function were 
used to obtain the predictions for the large distance asymptotic forms 

T / n \ 9RR T i n \ 9RG t i n \ 9GG 

where A := g rrQgg ~ Org- These are independent of the particle density. In contrast, for the 
true log-potential plasma system, charges +1 and +2, it is only the combination 

P+i,+i(^ 0) + Ap T +l +2 (x, 0) + 4pJ 2 +2 (z, 0) , 



r 2™2 



7r°ar 

relating to the charge-charge correlation which exhibits a density independent large x form. 
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A Appendix 



In [10] , Ishikawa, Okada, Tagawa and Zeng give the following Pfaffian identity for the Vander- 
monde determinant. 

Theorem A.l Suppose N is an even integer and let x = (x\,xi, . . . , xjv) be indeterminants. 
Then, 

■ (N/2 _ N/2\%-\ N 

= 11 fcn-Sm)- (66) 
m,n=l l<m<n<N 



Pf 



Ishikawa et al. prove this identity as a special case of a more general Pfaffian evaluation. A 
simpler proof, following that of the Vandermonde determinant identity, is outlined as follows: 

1. The left hand side of fl66]) is if x n = x m for any n ^ m, and thus the right hand side 
divides the left ( clS db polynomial in Q[xi,x 2 , ■ ■ ■ ,x N ]; 

2. The polynomials defined by the left and right hands sides of flBBT) are homogeneous and 
of the same degree, and thus their ratio is a non-zero rational number; 

3. This rational number is equal to 1 as seen by checking that the coefficients on the left 
and right hand side of any particular monomial are equal. 



This proof works more generally. If F(x) and G(x) are monic polynomials of degree N/2, then 
Pf 



(F(x n ) - F(x m )) (G(x n ) - G{x m j) 



I Xri. Xv 



m,n=l l<m<n<N 
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The following Theorem is similar in spirit to these Pfaffian evaluations, but is in a form 
which lends itself more naturally to calculations which arise in random matrix theory and 
related fields. 

Theorem A. 2 Suppose N is an even integer and let x = (x±,X2, ■ ■ ■ ,xn) be indeterminants. 
If 7To, 7Ti, . . . , 7Tjv-i are polynomials with deg7r n = n and leading coefficients ao, a±, . . . , ajv-i, 
then 

-l N 



Pf 



N/2-1 



N-l 



II ak n ( X n-Xm)- 

k=0 l<m<n<N 



- m,n=l 

In addition, we prove a confluent form of this identity. For each non-negative integer £ we 
define the differential operator 

(D is the identity operator). If L is a positive integer and 7r , tt±, . . . , ttln-i are polynomials 
with degp n = n and leading coefficients a ,a±, . . . , olat-i, we define the LN x LiV confluent 
Vandermonde matrix 



V 



DoKn-l( x r, 
Dx-Kn-AXr, 



_D L - 1 'K n - 1 {x m )_ 
The confluent Vandermonde identity has that 



m=l,2,...,7V; 
n = 1,2,... LN. 



(67) 



LN-l 



detV = Yl a k | ] (x„-x m ) L2 . 



fc=0 l<m<n<N 



Theorem A. 3 Suppose L is an integer, N is an even integer and 7To, 7Ti, . . . , 7Tln-i are poZ?/- 
nomials with degp n = n and leading coefficients a , a±, . . . , a^jv-i- Define, 



LN/2-1 



K(x,y)= 

^2n{x)7T2n+l{y) - ^2n+l {x)n 2n (y) ■ 



(68) 



n=0 



and /or eaca < £, m < L — 1 define 

LN/2-1 

D e K(x,y)D m = J2 D i iT2 n (x)D m n2 n+1 (y)-D i n2 n+1 (x)D m 7r2 n (y). 



n=0 



Define the antisymmetric LN x LiV matrix 



W = 



D°K{x m ,x n )D° D°K{x m ,x n )D l 
D 1 K(x m ,x n )D° D 1 K(x m ,x n )D 1 



D*K{x m ,x n )D L - x 
D x K{x m ,x n )D L - x 



D L - l K{x m ,x n )D° D L ~ l K{x m ,x n )D l ••• D L ~ x K{x m , x n )D 



L-l 



N 



m,n=l 
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Then, 



LN-l 



= n ^ n 



\L 2 



k=0 Km<n<N 



Clearly Theorem IA.2I follows from Theorem IA.3I by setting L = 1. 
Corollary A. 4 Let 



F(x,y) 



( y LN/2 _ x LN/2y. 

(y-x) 



and for each < £, m < L — 1 let 
D e F(x,y)D m 



1 



lim lim — — -- — F(w, z) 



(£ — l)!(m — 1)! \w^xz^y dw e dz 
and define the LN x LN antisymmetric matrix 



U 



D°F(x m ,x n )D° D^ixm.x^D 1 
D l F(x m) x n )D° D 1 F(x m , injD 1 



D°F(x m ,x n )D L ~ 1 
D 1 F{x m ,x n )D L - 1 



-l N 



D L ~ 1 F(x m , x n )D° D L - x F{x m ,x n )D x ■■■ D L ^ 1 F(x mi x n )D 



L-l 



J m,n=l 



Then, 



PfU 



n 



\X n X m j 



l<m<n<N 

Proof of Theorem IA.3I We define the LN x LN matrix 



J 



1 
-1 



1 
-1 



1 
-1 



It is an easy exercise to check that Pf J = 1. A similarly easy calculation shows that W = VJV T , 
and thus, using properties of the Pfafnan, PfW = det VPfJ = det V. □ 
Proof of Corollary IA.41 Suppose we permute and relabel n , tti, . . . , kln-u so th & t it is not 
necessarily the case the deg7r n — n. If a is the permutation of 0, 1, ... , LN — 1 which we used 
to permute our polynomials, then it is easy to see that V, formed as in fl67|) . but using our 
reordered and relabeled polynomials will satisfy 



LN—l 



det V = sgncx jQ ] [ 



{x n Xiri) , 



fc=0 Km<n<N 
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and if we define K(x,y) as in (168)) . but with our reordered and relabeled polynomials and 
defining W with this 'new' K, we have 



pfw = sgna n ^ n 



LN—l 

•En •^rn ) 



k=0 l<m<n<N 

Now, 

(„.LN/2 _ ^,LN/2\2 LN/2-1 
F ( x ,y) = ^—, ^_J_= y y LN-l-n x n_ y n x LN-l-n 

(y-x) ^ 

Thus, if we define 

K2n+i{x) = x LN ~ l ~ n and 7r 2n {x) = x n ; n = 0, 1, . . . , LN/2 - 1, 

we have 

LN/2-1 

F(x,y) = ^ n 2n(x)7r 2n+1 (y) - vr 2n+ i(x)7r 2 „(y) = K(x,y). 

71=0 

Finally, the signature of the permutation which orders ir , 7Ti, . . . , kln-i by degree can be shown 
to always be 1 (since the reverse permutation on LN/2 elements has the same signature as the 
perfect shuffle on LN elements). □ 
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